Ocean circulation contributes to genetic connectivity of limpet populations at deep‐sea hydrothermal vents in a back‐arc basin

Abstract For endemic benthos inhabiting hydrothermal vent fields, larval recruitment is critical for population maintenance and colonization via migration among separated sites. The vent‐endemic limpet, Lepetodrilus nux, is abundant at deep‐sea hydrothermal vents in the Okinawa Trough, a back‐arc basin in the northwestern Pacific; nonetheless, it is endangered due to deep‐sea mining. This species is associated with many other vent species and is an important successor in these vent ecosystems. However, limpet genetic diversity and connectivity among local populations have not yet been examined. We conducted a population genetics study of L. nux at five hydrothermal vent fields (maximum geographic distance, ~545 km; depths ~700 m to ~1650 m) using 14 polymorphic microsatellite loci previously developed. Genetic diversity has been maintained among these populations. Meanwhile, fine population genetic structure was detected between distant populations, even within this back‐arc basin, reflecting geographic distances between vent fields. There was a significant, positive correlation between genetic differentiation and geographic distance, but no correlation with depth. Contrary to dispersal patterns predicted by an ocean circulation model, genetic migration is not necessarily unidirectional, based on relative migration rates. While ocean circulation contributes to dispersal of L. nux among vent fields in the Okinawa Trough, genetic connectivity may be maintained by complex, bidirectional dispersal processes over multiple generations.

vent species are benthic; therefore, dispersal that connects the populations and maintains genetic diversity of metapopulation is mediated by pelagic larvae or mobile juveniles (Vrijenhoek, 2010).For benthos endemic to hydrothermal vent environments, larval recruitment is critical for population maintenance and colonization via migration among populations, as vent fields are patchy on the seafloor (Adams et al., 2012).Long-distance planktonic larval dispersal contributes to genetic diversity and connectivity among vent-field populations, and may help to mitigate degradation caused by disturbances.Most larvae die due to predation and lethal stresses during planktonic dispersal.
Surviving larvae that manage to reach other hydrothermal vent fields, settle and grow.Larvae mature and reproduce to leave a genetically traceable footprint in the local population.Long-distance larval dispersal likely occurs stochastically in stepping-stone fashion over many generations (Geiger & Thacker, 2005;Vrijenhoek, 1997).Therefore, recruitment among hydrothermal vent fields appears to have contributed to expansion of species distributions at a historical timescale, and maintenance of metapopulations at an ecological timescale.Mitarai et al. (2016) predicted the probability of larval dispersal among hydrothermal vents using an ocean circulation model, and suggested possible connections between hydrothermal vent fields.
Lecithotrophic larvae, which are dependent on nutrition derived from egg yolk, appear to have shorter planktonic larval durations and generally more limited dispersal potential than planktotrophic larvae, which are able to feed while dispersing (Lutz et al., 1984;Lutz & Kennish, 1993;Vrijenhoek, 2010).Nonetheless, larval nutrition type alone may not determine actual larval dispersal distance, as larval duration in lecithotrophic species may be prolonged by production of large eggs or arrested development (Vrijenhoek, 2010).Moreover, empirical assumptions about larval dispersal may not apply to deep-sea species due to colder water temperatures.Lower temperatures and slower metabolism can prolong larval duration, although they also increase mortality due to increased predation and severe environmental conditions (O'Connor et al., 2007).Most hydrothermal vent gastropods appear to have lecithotrophic larvae (Tyler et al., 2008).It is likely that Lepetodrilus is no exception (Craddock et al., 1997;Johnson et al., 2006Johnson et al., , 2008;;Plouviez et al., 2009;Vrijenhoek, 2010), though the egg size is suggestive of planktotrophic larval development (Tyler et al., 2008).
Nevertheless, a low rate of successful larval dispersal appears sufficient to maintain connectivity among sites (Lowe & Allendorf, 2010).
Population genetic studies using genetic markers are essential for estimation of population connectivity among geographically isolated sites.Since, genetic connectivity is influenced by larval dispersal because benthic animals disperse during the larval stage, vent gastropod species are genetically differentiated by geographic distance between regions, that is, horizontal distance (Ifremeria nautilei; Tran Lu et al., 2022) or between depths, that is, vertical distance (Lepetodrilus concentricus, formerly Lepetodrilus sp.; Roterman et al., 2016).However, there is little genetic differentiation within regions in Lepetodrilus nux (Nakamura et al., 2014).In the case of Shinkailepas, a planktotrophic limpet, there is no significant differentiation between regions separated by 1000 km (Yahagi et al., 2017(Yahagi et al., , 2019(Yahagi et al., , 2020)).In vent gastropods, genetic connectivity has been studied using one to several genetic loci from mitochondrial and/or nuclear gene(s) (Breusing et al., 2023;Nakamura et al., 2014;Thaler et al., 2011;Yahagi et al., 2017Yahagi et al., , 2019Yahagi et al., , 2020)).Meanwhile, high-resolution genetic markers for vent species are expected to detect fine genetic connectivity, and some studies assessed connectivity using microsatellites or single-nucleotide polymorphisms (SNPs) (Roterman et al., 2016;Tran Lu et al., 2022).Studies of regional genetic connectivity based on multilocus genotypes between hydrothermal vents are scarce, due to the difficulty of collecting samples.Nonetheless, genetic markers have helped to illuminate population structure and dynamics of larval dispersal in population maintenance of deep-sea hydrothermal vent species (Roterman et al., 2016;Thaler et al., 2011).In addition, next-generation sequencing has facilitated development of polymorphic microsatellite markers for population genetic analysis.
The Okinawa Trough is a back-arc basin along the Ryukyu Archipelago in the northwestern Pacific Ocean.This trough is approximately 1000 km in length and more than 10 hydrothermal vent fields have been discovered in it to date (Fujikura et al., 2008;Miyazaki et al., 2017;Nakamura et al., 2015;Watanabe et al., 2010).
However, some of these hydrothermal vent fields have been disturbed by scientific drilling or mineral exploration (Kawagucci et al., 2013).Populations of hydrothermal, vent-endemic lepetodrilid limpets, Lepetodrilus nux, maintain considerable genetic diversity and are well-mixed among vent fields in the Okinawa Trough, as illustrated using the mitochondrial COI gene (Nakamura et al., 2014).In this study, we documented genetic diversity and connectivity of L. nux using 14 polymorphic microsatellite loci developed using nextgeneration sequencing (Nakajima et al., 2018).Now it is possible to gather more detailed genetic diversity and connectivity data using microsatellites, to better understand mechanisms of population maintenance in L. nux.Lepetodrilus nux is an endangered species (Molloy et al., 2020), but widely distributed at hydrothermal vents in the Okinawa Trough (Nakamura et al., 2014;Okutani et al., 1993;Sasaki et al., 2003).The genus Lepetodrilus is one of most diversified gastropod groups at deep-sea hydrothermal vents (Desbruyères et al., 2006), and some Lepetodrilus species also colonize cold seeps, wood falls, and whale carcasses (Johnson et al., 2008).To consider the timescale for connectivity, we compared genetic connectivity to potential larval dispersal from the ocean circulation model (Mitarai et al., 2016).This model simulated dispersal times for hydrothermal vent sites using 55 years of data.Therefore, the model is appropriate for quantitative comparisons of population genetic connectivity and current ocean circulation patterns.

| Collection of samples and scoring of microsatellite genotypes
Specimens of L. nux were collected at five hydrothermal vent fields in the Okinawa Trough (maximum geographic distance: ~545 km, depth range: ~700 m to ~1650 m) (Table 1), using the remotely operated vehicle, Hyper-Dolphin, during cruises NT11-20 (September 29, 2011-October 12, 2011) and NT13-22 (November 7-19, 2013) of the R/V Natsushima, conducted by the Japan Agency for Marine-Earth Science and Technology.All specimens were preserved in 99.5% ethanol and then transferred to an onshore laboratory.Genomic DNA was extracted using a DNeasy Blood & Tissue Kit (Qiagen).

| Documenting genetic diversity
Numbers of alleles, values of observed and expected heterozygosity (H O and H E , respectively), and a deviation index (F IS ) from HWE were also evaluated with GenAlEx ver.6.501 (Peakall & Smouse, 2006).
For estimation of genetic diversity, we calculated allelic richness at each vent field, normalized for differences of sample size using FSTAT ver.2.9.3.2 (Goudet, 1995).The expected H E value was used as an index of genetic diversity.The probability of recent genetic bottlenecks was examined using BOTTLENECK ver.1.2.02 (Piry et al., 1999) based on assumptions for both the infinite allele model (IAM) and the two-phase model (TPM; 30% IAM and 70% stepwise mutation model).For this analysis, we adopted Wilcoxon's signedrank test with 1000 replications.The probability of identify (p ID ) was also calculated with these 14 microsatellite loci, using GenAlEx to estimate the resolution of individual identification.

| Genetic clustering approaches
Population genetic structure was inferred using STRUCTURE ver.2.3.4 (Pritchard et al., 2000).STRUCTURE analysis implements a Bayesian clustering algorithm to assign genotypes to clusters that minimize HWE and linkage disequilibrium.Ten replicate runs were conducted for each K between 1 and 7 with location prior information using the admixture model and assuming correlated allele frequencies (Falush et al., 2003).
Each run consisted of 1,000,000 Markov chain Monte Carlo (MCMC) replications after burn-in with 100,000 iterations.Optimal K was determined using the method of Evanno et al. (2005), as implemented in STRUCTURE HARVESTER (Earl & vonHoldt, 2012).Run data were merged with CLUMPAK (Kopelman et al., 2015).Discriminant Analysis of Principal Components (DAPC) was conducted in R ver.3.4.4(R Core Team, 2018) using the package, adegenet ver.1.3-9.2(Jombart et al., 2010) to represent genetic patterns of each individual.This clustering method does not assume a particular model; therefore, it is free of assumptions about HWE and linkage disequilibrium.Clusters were predefined for each site and different distances among clusters were demonstrated in a scatterplot of individuals.

| Directional relative migration rates between vent fields
We employed the divMigrate function (Sundqvist et al., 2016) in diveRsity ver.1.9.90 (Keenan et al., 2013) in R to estimate relative migration rates between vent fields, with the classical measurement method of genetic differentiation using allele frequency data based on the infinite island model.Statistical methods used Jost's D (Jost, 2008).The maximum migration rate between sites is normalized as 1, and other pairwise sites are reported as relative migration rates (<1).An asymmetrical test to estimate unidirectional migration was conducted and 95% confidence intervals were calculated from 15,000 bootstrap samples.In addition, we tested the significance of the directionality difference for both horizontal (geographic distance) and vertical migration (difference of depth), using values of relative migration rates.

| Estimation of potential long-distance dispersal based on an ocean circulation model
Estimation of successful dispersal among five vent fields was based on the ocean circulation model from the Regional Ocean Modeling System (ROMS) (Mitarai et al., 2016), which covers the five vent fields in this study.Planktonic larval duration (PLD) was set to the water temperature for each depth in western Pacific vent fields, estimated by Mitarai et al. (2016) following the previous mean PLD data of 69 marine species, including planktotrophic and lecithotrophic larval species (O'Connor et al., 2007).This is because larval metabolism can be accelerated by higher water temperature, so larvae have shorter PLDs closer to the surface.For example, larval dispersal periods of six species of lecithotrophic mollusks were estimated as 15.4 days at 10°C and 8.6 days at 20°C.In contrast, those of planktotrophic mollusks were estimated for 10 species at 50.2 days and 27.3 days, respectively (O'Connor et al., 2007) (Figure S2).The ocean circulation model shows the number of successful dispersal events of passive larvae in 100 million independent dispersal events from each vent site (Mitarai et al., 2016).We calculated the probability based on the number of successful dispersal events that can connect a focal population to its natal population per 100 million dispersal events.The probability was assessed per generation and per two successive generations, to simplify connectivity among vent fields.The probability of two successive generations was assumed to represent steppingstone connectivity among vent fields.The probability was calculated for depths of 100, 300, 500, 700, and 1000 m.The probability of recruitment at MEK was not considered at 1000 m, because the vent occurs at ~700 m (Table 1).Furthermore, the probability of dispersal was estimated in both directions, from northeast to southeast and from southeast to northeast, as relative migration rates.

| Genotyping of Lepetodrilus nux
Using 14 microsatellite loci (Table S1), we scored multilocus microsatellite genotypes of 167 individuals collected from five hydrothermal vent fields in the Okinawa Trough (Table 1).Complete multilocus genotypes were obtained for 158 individuals for use in subsequent analyses.The probability of identity for each vent field ranged from 2.3 × 10 −17 to 1.5 × 10 −15 .There were no replicated multilocus genotypes; therefore, all individuals in our sample possess unique multilocus genotypes, confirming the utility of these 14 loci.

| Genetic diversity and bottlenecks in populations
Mean allelic richness (A R ), normalized for differences in sample size across 14 loci per vent field, ranged from 8.94 to 9.93 (Table 1).
Sixty private alleles (P A ) were detected in 14 loci at five vent fields (Table S2).Mean observed and expected heterozygosities (H O and H E ) ranged from 0.674 to 0.751 and from 0.731 to 0.778 (Table 1), respectively.The mean deviation index (F IS ) from HWE ranged from −0.022 to 0.117 (Table S2).Although three vent fields (IZH, IRK, and HTK) appear to have experienced a recent bottleneck under the infinite allele model (IAM) (p < 0.05), no populations displayed such a bottleneck under the two-phase model (TPM) (p < 0.05) (Table 1).

| Genetic differentiation among hydrothermal vent fields
Hierarchical AMOVA indicated significant genetic differentiation among vent fields (F ST = 0.009, p = 0.001), and pairwise F ST and G ′′ ST values ranged from −0.003 to 0.023 and from −0.016 to 0.098, respectively (Table S3).Significant genetic differentiation of both F ST and G ′′ ST was confirmed between distant locations, whereas the closest hydrothermal vents did not show significant differentiation (Figure 1a).Supporting the results of F ST and G ′′ ST , STRUCTURE analysis with the LOCPRIOR model showed differences in the Q-matrix among sites, especially those most separated.The Evanno method indicated that the most probable number of genetic populations was three, as the highest indicated ΔK was K = 3 (Figure S1a).However, genetic structure is not strong because each cluster is largely shared between vent fields (Figure 1a and Figure S1b).DAPC also documented differences in cluster composition between vent fields.In that analysis, the first and second discriminant functions explained 51.56% and 18.95% of the variance, respectively.Clusters, IRK and HTK, are especially isolated from the other three overlapping population clusters (Figure 1b).Those three clusters are closest to the IRK cluster, and the IRK cluster is most similar to the HTK cluster.
There is no obvious differentiation in relation to depth (from ~700 m to ~1650 m) in this species.Isolation-by-distance was significant (p = 0.010 in F ST , p = 0.017 in G ′′ ST ) (Figure 2a,b), but isolation-by-depth was not confirmed (p = 0.350 in F ST , p = 0.336 in G ′′ ST ) (Figure 2c,d).S3.INF, and IZH exceeded 0.5, regardless of direction (Figure 3b).In addition, the migration rate was higher from northeast to southwest than the rate from southwest to northeast (Figure 3c); however, the two migration rates are not significantly different.Furthermore, mutual migration patterns between deep and shallow vent fields were detected (Figure 3d).

| Comparison of isolation-by-distance with other species
For vent-only populations, we detected isolation-by-distance in G.
platifrons, analyzed by genome-wide SNPs (Figure S3a) regardless of planktotrophic nutrient larval type and with longer PLD.In contrast, the other three species (P.subglabra, B. nipponica, and S. crosnieri) and G. platifrons analyzed by mtDNA did not show isolation-by-distance (Figure S3a).However, B. nipponica showed relatively high values of F ST (ranging from 0.0075 to 0.3431), whereas S. crosnieri showed low F ST values ranging from −0.0187 to −0.0040.
For both vent and seep populations, we did not detect significant isolation-by-distance in any of the four vent-related species.In terms of SNPs, G. platifrons also did not show isolation-by-distance due to inclusion of seep populations, including the Jiaolong Ridge population in the South China Sea, which is relatively near the Okinawa Trough (Figure S3b).Shinkaia crosnieri did not exhibit isolation-bydistance, but it displayed large genetic differentiation between vent and seep populations with higher F ST values ranging from −0.0187 to 0.2301 (Figure S3b).

| Probability of dispersal based on the ocean circulation model
The probability of successful dispersal is relatively higher in deepwater (Figure 4).Maximum probabilities are 0.00568% at 700 m and 0.00860% at 1000 m per generation, assuming 100 million independent events.In addition, maximal probabilities of two-generation dispersal were also higher in deeper water, 0.00502% at 700 m and 0.01312% at 1000 m.However, the probability of two-generation dispersal in shallow water showed a lower value (0.00051%) at 100 m depth.Northward dispersal was relatively higher regardless of depth, both per generation and per two successive generations.Meanwhile, the relative migration rate using a genetic approach shows nonsignificant bidirectional migration patterns (p = 0.182) (Figure 3c).

| DISCUSS ION
Lepetodrilus nux populations in the Okinawa Trough show that there is no recent bottleneck effect at any vent field in this study.Genetic There is no significant difference between migration from southwest to northeast and from northeast to southwest.In addition, gene flow occurs regardless of depth (no isolation-by-depth).Although larvae appear to have buoyancy and vertical swimming ability, active behavior may not be sufficient to migrate long distances (over a few hundred kilometers).This results in limited larval dispersal range (Craddock et al., 1997;Lutz et al., 1984).Population genetic data and the ocean circulation model predict that bidirectional migration patterns are maintained by deeper currents over multiple generations, and the migration per generation mainly occurs between geographically proximal sites.These patterns are complex, due to eddies and countercurrents around the Kuroshio Current (Nakamura et al., 2013).Larvae of gastropods, such as limpets, tend to remain near the bottom or to wander with bottom currents (Mullineaux et al., 2005(Mullineaux et al., , 2013)).Roterman et al. (2016)  Megaplumes associated with frequent seafloor volcanic eruptions promote larval transport and such events may maintain genetic connectivity among sites, compensating for differences in depth (Matabos et al., 2008).Therefore, some larvae may be transported by upper currents; however, successful dispersal must be rare.As a counteracting effect, PLD becomes shorter and limits dispersal distance closer to the surface than at greater depth due to accelerated larval metabolism caused by increasing water temperature, compensating for increasing current speed (Mitarai et al., 2016).These new insights may be helpful to establish conservation strategies for hydrothermal vent gastropods against drilling for sea floor resources (Nakamura et al., 2018).
We conclude that connectivity of L. nux in the Okinawa Trough is limited by geographic distance (within ~545 km), independent of differences in depth.Previous studies have also reported limited genetic connectivity and differentiation of lecithotrophic vent gastropods.Restricted contemporary gene flow was found in Ifremeria nautilei between the Manus Basin and the Fiji-Lau area (Kojima et al., 2000;Thaler et al., 2011;Tran Lu et al., 2022).Moreover, there is no definite genetic break between Fiji and Lau and due to unknown processes, genetic connectivity does not reflect geography (Thaler et al., 2011).Plouviez et al. (2019) further reported that Lepetodrilus The probability of dispersal between hydrothermal vent fields estimated using the ocean circulation model (Mitarai et al., 2016).The probability of larval transport between hydrothermal vent populations per generation and per two successive generations at each depth (100, 300, 500, 700, and 1000 m).Data related to MEK (~700 m depth) could not be determined at 1000 m.Furthermore, dispersal probability is estimated in each direction.The direction between IZH and INF was considered a case of relative migration rates.Plot at 0 km of geographic distance means probability of self-recruitment or return at each vent field.The polynomial regression curve was produced using ggplot2 package in R.
aff. schrolli showed no genetic structure within basins, using 42 nuclear DNA markers, but that it showed large genetic differences between the Manus and Lau Basins.Vent-related mussels of the genus Bathymodiolus, a planktotrophic larval taxon that may be able to remain in the water column or on the surface for up to 1 year (Arellano et al., 2014;Arellano & Young, 2009), showed high genetic similarity among regions (Breusing et al., 2015).The planktotrophic red blood limpet, Shinkailepas myojinensis, shows no genetic differentiation between the Okinawa Trough and the Izu-Ogasawara Arc (~1350 km) based on an analysis of the mitochondrial COI gene, as its larvae are expected to have an extraordinarily long planktotrophic larval period of more than a year, and surface temperature limits their distribution (Yahagi et al., 2017).
A population genomic analysis of Gigantidas platifrons, which is taxonomically close to Bathymodiolus, also showed no significant genetic differentiation among populations in the Okinawa Trough, based on an analysis of genome-wide SNPs (Xu et al., 2018).Larvae of G. platifrons are planktotrophic, with longer larval duration than lecithotrophic larvae (Xu et al., 2018).Nevertheless, isolation-bydistance was detected in vent populations of G. platifrons analyzed by genome-wide SNPs.Low genetic differentiation and bidirectional migration patterns were also shown in B. nipponica and S.
crosnieri in the Okinawa Trough, analyzed by genome-wide SNPs (Xu et al., 2021(Xu et al., , 2024)), and there is no isolation-by-distance in these two species (Figure S3a).In addition, seep populations contribute to isolation of vent populations due to heterogeneity that disturbs  (Vrijenhoek, 2010).This larval trait largely determines larval dispersal capacity and genetic differentiation among sites.Although the genus Lepetodrilus is considered lecithotrophic (Vrijenhoek, 2010), larvae of Lepetodrilus are potentially planktotrophic, judging from egg size and fecundity (Tyler et al., 2008) (Minami & Ohara, 2017), located between IZH and IRK, may offer genetic connectivity in the Okinawa Trough.Further study of population genetics with more populations may enable us to better understand connectivity and migration patterns of this species.
The ocean circulation model by ROMS estimated that successful migration between hydrothermal vents tends to occur at depths of 700-1000 m.Surface currents may strongly disperse larvae, reducing the rate of successful recruitment.Ocean circulation at greater depths may contribute to short-distance dispersal to neighboring vent fields at this geographic scale.In one generation, the probability of larval transport is indeed higher at shallower depths than at 1000 m over longer distances, for example, over 400 km, whereas in two generations, it becomes greater at 700 and 1000 m depths, than in shallower water at all geographic distances (Figure 4).have been reported between genetic approaches and the ocean circulation model (Breusing et al., 2023;Moody et al., 2019).The timescale difference between genetic approaches via multiple generations with evolutionary history and the ocean circulation model may create such discrepancies.Temporal variation in larval supply potentially has important consequences for population genetic structure (Mullineaux et al., 2005).Historic colonization and larval physiology are crucial for genetic connectivity and these factors vary even among related species (Breusing et al., 2023).However, the L. nux PLD is unknown and should be determined to better define the actual dispersal range of this species.PLD variation may contribute to recruitment at both natal and distant vent fields.
Further studies of larval behavior, actual dispersal depth, and larval physiology will allow realistic estimates of larval dispersal and survivorship.In addition, "ghost populations" with unknown effec- . Genetic differentiation indices between vent fields, pairwise F ST , and G ′′ ST , were calculated using GenAlEx.The significance of each F ST value was tested with 999 permutations.Furthermore, the G-statistic method was used to calculate pairwise G ′′ ST , which is effective for estimating genetic differentiation with a small number of populations (Meirmans & TA B L E 1 Location information and genetic indices for Lepetodrilus nux at each sampled vent field in the Okinawa Trough. is the number of samples analyzed.A R is allelic richness, the standardized index of genetic diversity.H O and H E are the mean observed and expected heterozygosities, respectively.F IS is the fixation index (inbreeding coefficient), deviation from HWE. IAM and TPM are significant values of the infinite allele model and two-phase model in the analysis using BOTTLENECK, respectively.p ID is the probability of identify for each vent field.Hedrick, 2011), and the significance of each value was tested with 999 permutations.Isolation-by-distance was examined using the Mantel test in GenAlEx to determine whether the correlation between pairwise F ST /(1 -F ST ) and G ′′ ST /(1 -G ′′ ST ) values and direct geographic distances (km) between vent fields was significant.Isolation-by-depth was also analyzed by the same method.Furthermore, the isolationby-distance pattern was compared to patterns of other vent species, including vent-related species inhabiting both vent and seep sites, in the Okinawa Trough based on data published previously.Isolationby-distance in published data was calculated from the genetic differentiation based on pairwise F ST or Φ ST , and geographic distance was calculated from coordinates of vent sites.Published data are derived from analyses of the abyssochrysoid gastropod, Provanna subglabra divMigrate analysis revealed gene flow among vent fields.A pattern of isolation-by-distance was detected based upon F ST and G ′′ ST values.The highest migration rate occurred between INF and IZH, and the lowest rate was from HTK to MEK.These results also reflect geographic distance (Figure 3a).All migration rates among MEK, F I G U R E 1 (a) Genetic clusters among five vent fields based upon the Q-matrix, estimated by STRUCTURE analysis.The number of clusters are three (K = 3), and the optimal K was determined using the ∆K method of Evanno et al. (2005) (Figure S1a).Colors of the pie charts show each genetic cluster.Bar plots for each individual are shown in Figure S1b.Arrows show non-significant differentiation in both F ST and G ′′ ST between vent fields.Significance levels and pairwise F ST and G ′′ ST values are shown in Table S3.Contour lines indicate sea depths at intervals of 500 m.(b) Discriminant analysis of principal components (DAPC) shows scatterplots with prior site information.Each plot explains multilocus genotypes analyzed using 14 microsatellites.F I G U R E 2 Isolation-by-distance (a, b) and isolation-by-depth (c, d) result from the F ST /(1 -F ST ) (a, c) and G ′′ ST /(1 -G ′′ ST ) (b, d) values of Lepetodrilus nux in the Okinawa Trough.Pairwise F ST and G ′′ ST values are shown in Table

F
I G U R E 3 (a) Relative migration rates among populations of Lepetodrilus nux at five vent fields in the Okinawa Trough estimated with divMigrate.The most plausible migration rate was assumed as 1 (from INF to IZH), and relative migration patterns were output for other vent fields in pairwise fashion.An asymmetrical test revealed that all unidirectional migration patterns were non-significant at a 95% threshold.(b) Relative migration rates between vent fields were over 0.5.Relative migration rates in each direction by geographic distance (c), and by depth (d), respectively.Migration from IZH to INF was counted as southwest to northeast in this analysis, and vice versa.
differentiation and clustering based on microsatellites reflect geographic distance, in contrast to results using the mitochondrial COI gene.Nevertheless, these genetic data appear to support a low level of genetic differentiation of L. nux in the Okinawa Trough.Genetic connectivity has been maintained, even at MEK, which is relatively shallow, compared to the other four vent fields.Geographic distance tends to prevent vent-to-vent larval transport, but there is no gradual genetic differentiation by depth.Relative migration analysis showed that the highest migration rate occurs between INF and IZH.Three vent fields (MEK, INF, and IZH), which are geographically close, show bidirectional migration.
reported that L. concentricus (formerly Lepetodrilus sp.) populations at the East Scotia Ridge (ESR) are significantly differentiated from a population in the South Sandwich Island Arc (SSI), only ~95 km from ESR, but there is a large difference in depth between ESR and SSI (2396-2645 m vs. 1434 m depth).Large depth differences may promote genetic differentiation, or deeper environments may result in adaptive isolation of Lepetodrilus populations.Nevertheless, differences in depth are not a barrier to genetic connectivity of L. nux in the Okinawa Trough.
Isolation-by-depth is non-significant; therefore, stochastic recruitment occurs from ~700 m (MEK) to ~1650 m (IRK) depth.However, the Kuroshio Current is very strong along the Okinawa Trough and the ocean circulation model tends to support large potential migration from southwest to northeast.Contrary to expectation, migration from northeast to southwest appears to occur in adjacent areas.Ocean circulation via the Kuroshio Current from southwest to northeast appears likely as a driving force of long-distance dispersal.Genetic migration rates show bidirectional migration without significant directionality.Discrepancies tive population sizes derived from unsampled or unknown hydrothermal vent fields may support migration, contributing to genetic connectivity.The existence of Lepetodrilus species in other regions in the northwestern Pacific may reveal previously unknown populations.Moreover, investigation of hydrothermal vent-endemic invertebrates using genotypic big data, such as genome-wide, single-nucleotide polymorphisms with quantitative trait loci may clarify the distribution and adaptation of this species/genus in the northwestern Pacific.